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Abstract. A contact interaction analysis is presented to search for new phenomena 
beyond the Standard Model in deep inelastic e^p — > hadrons scattering. The 
data are collected with the HI detector at HERA and correspond to integrated 
luminosities of 0.909 pb~^ and 2.947 pb~^ for electron and positron beams, respec- 
tively. The differential cross sections da/dQ"^ are measured in the range between 
160 GeV^ and 20,000 GeV^. The absence of any significant deviation from the 
Standard Model prediction is used to constrain the couplings and masses of new 
leptoquarks and to set limits on electron-quark compositeness scales and on the 
radius of light quarks. 
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1 Introduction 



The e^p collider HERA, in which 27.5 GeV (26.7 GeV in 1993) leptons colhde with 820 GeV 
protons, provides access to an as yet unexplored mass domain for the discovery of new particles. 
Several extensions of the Standard Model postulate either new fermions or new bosons, such 
as leptoquarks and additional gauge bosons. Signals for physics beyond the Standard Model 
can be discovered either directly or indirectly. 

Common to all direct searches is the s channel formation of a new heavy resonance at a mass 
Mx = \/x~s, which is, however, limited by the available centre of mass energy of a/s ~ 300 GeV. 
The scaling variable x is the momentum fraction of the proton carried by the struck quark. 
Results of direct searches for new heavy bosons at HERA have been published recently by the 
experiments HI and ZEUS §. 

The search for new bosons or eg compositeness can be considerably extended beyond the 
kinematic production limit through the study of indirect effects from virtual particle exchange. 
Such effects may become observable as deviations from the Standard Model expectation at high 
momentum transfers Q^. 

The analysis presented here combines all data collected with the HI experiment at HERA 
during 1993 and 1994 with electron and positron beams. They represent an increase in inte- 
grated luminosity by a factor of ~ 9 compared to the previous data sample 



2 Phenomenology of (e e) {q q) Contact Interactions 

New currents or bosons may produce indirect effects through the interference of a wzrtwa/ particle 
exchange with the 7 and Z fields of the Standard Model. For particle masses well above the 
available production energy, such indirect signatures may be investigated by adding general 
contact interaction terms to the Standard Model Lagrangian. Sufficiently heavy particles X 
cease to propagate and thus new contact terms and modified vertices arise from 'contracting' the 
particle propagators to an effective 4-fermion point-like interaction. The separate dependence 
of s, t and u channel amplitudes on couplings gx^i, / to states i, / and mass Mx reduce to the 
dependence on effective couplings with dimension [mass~^] 

"^'^ - Ml • 

The most general chiral invariant neutral current contact interaction Lagrangian can be 
written in the form |0] 

q = u,d 

+ Vrl {eRlf,eR){qLl^qL) + r]1^R {eRjf,eR){qRYqR)} , 

where the indices L and R denote the left-handed and right-handed fermion helicities and the 
sum extends over up and down quarks and antiquarks q. 

Although contact interactions have been originally proposed in the context of composite 
leptons and quarks ^, this ansatz can be easily applied to other new phenomena by an 
appropriate choice of the coupling coefficients rjif. 

Leptoquarks are colour triplet bosons of spin or 1, carrying lepton (L) and baryon (B) 
number and fractional electric charge. They couple to lepton-quark pairs and appear in almost 



5 



all extensions of the Standard Model which try to establish a connection between leptons and 
quarks @, |^. Leptons and quarks may be either arranged in common multiplets, like in Grand 
Unified Theories or superstring motivated Eq models, or they may have a common substructure 
as in composite models. A fermion number F = L + 3 B is defined, which takes the values F = 2 
for leptoquarks coupling to e~ q and F = for leptoquarks coupling to e" q. For positrons the 
fermion number F changes by two units. Consequently F = 2 leptoquarks are easier accessible 
in e~ p scattering, while positron beams are more sensitive to F = leptoquarks, since at 
moderate x values quarks are more abundant in the proton than antiquarks. 

The notation, the contact interaction coefficients rjif and the fermion number assignment 
for leptoquarks with mass Mlq and coupling A are given in Table |^ (from ref. 0). The only 
unknown is the ratio M^q/X. Note that vector leptoquarks have positive coupling coefficients, 
while scalar leptoquarks have negative coupling coefficients, being a factor of 2 smaller in 
magnitude. 

In the Standard Model the fundamental particles - leptons, quarks and gauge bosons - are 
assumed to be pointlike. A possible fermion compositeness or substructure can be expressed 
through rjif = ±(7^/A^^, where the signs indicate positive and negative interference with the 
Standard Model currents, g is the coupling strength conventionally choosen as g'^/4:7T = 1 and 
A is the compositeness scale. 

3 The HI Detector 

A detailed description of the HI detector can be found elsewhere [Q. Here only those compo- 
nents are described, which are relevant for the present analysis. 

The lepton energy and angle is measured in a finely segmented liquid argon (LAr) sampling 
calorimeter covering the polar angleQ range 4° < 6 < 153° and all azimuthal angles. It consists 
of a lead/argon electromagnetic section with a thickness varying between 20 and 30 radiation 
lengths and a stainless steel/argon section for the measurement of hadronic energy fiow, which 
offers in total a containment varying from 4.5 up to 8 interaction lengths. Electron energies are 
measured with a resolution of a{E)/E ~ 12 %/\fE © 1% and hadron energies with a{E)/E ~ 
50 %/\/E(B2%. The absolute energy scales are known to 3% and 5% for electrons and hadrons, 
respectively. The angular resolution of the scattered lepton measured from the electromagnetic 
shower in the calorimeter is ~ 7 mrad. A lead/scintillator electromagnetic backward calorimeter 
extends the coverage at large angles (155° < 9 < 176°). The instrumented iron flux return yoke 
is used to measure the leakage of hadronic showers. 

Located inside the calorimeters is a tracking system, which consists of central drift and 
proportional chambers (25° < 9 < 155°), a forward track detector (7° < 6* < 25°) and back- 
ward proportional chambers (155° < 9 < 175°). The tracking chambers and calorimeters are 
surrounded by a superconducting solenoid coil providing a uniform field of 1.15 T within the 
tracking volume. 

The luminosity is determined from the rate of the Bethe-Heitler process ep — >■ epj mea- 
sured in a luminosity monitor far downstream the electron direction. The systematic errors 
of the integrated luminosity vary between 1.8% (1994 data) and 4.5% (1993 data). 



^The incoming proton moves in the +z direction with polar angle 6 — 0°. 
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4 Data Selection and Analysis 



The analysis is based on a purely inclusive measurement of the final state lepton in deep 
inelastic neutral current events p hadrons. All quantities of event properties are 

determined from the calorimeters alone, tracking information is only used to get the primary 
vertex position. The kinematic variables Q^, the negative squared momentum transfer, and the 
scaling variable y are derived from the scattered lepton energy E'^ and polar angle 6e 

= 4E,e:cos2|, 
y = l-^sin^^, 

where is the lepton beam energy. These quantities are related to the Bjorken scaling variable 

X hj Q"^ = xy s. 

The data were taken in three periods: (i) during 1993 with 26.7 GeV electrons, (ii) during 
1994 with 27.5 GeV electrons, and (iii) during 1994 with 27.5 GeV positrons. The event 
selection is similar to the previous leptoquark analysis [0 with the following requirements: 

1. The transverse energy of the scattered lepton has to exceed -^i^e > 8 GeV. 

2. The polar angle of the scattered lepton has to be within the LAr calorimeter acceptance 
10° < 0e < 150°. 

3. A primary vertex has to be reconstructed within | z^ertex ~ i^) \ < 35 cm of the nominal 
interaction point (z). 

4. The energy-longitudinal momentum conservation \ J2 [E — Pz) — 2 Ee \ < 10 GeV must 
be satisfied, where the sum extends over all detected particles or energy clusters. 

5. The event has to be balanced in transverse momentum | p™^ \ < 15 GeV. 

6. The scaling variable y has to fullfil y < 0.8. 

Requirements (1) and (2) assure that the kinematic quantities are well measured in the LAr 
calorimeter. is always well measured, while the resolution in y degrades at low values of 
y < 0.1. The requirements (3) and (5) suppress beam-wall and beam-gas background. The 
criteria (4) and (5) provide a good containment of the final state particles and reject events 
with a hard photon radiated from the initial state lepton and photoproduction events with a 
misidentified lepton in the LAr calorimeter. For events with all final state particles detected 
(except the proton remnants) one expects [E — Pz) ~ 2 Eg. Requirement (6) is introduced to 
avoid the region affected by large radiative corrections and to suppress photoproduction events 
with a misidentified lepton. The remaining contamination from photoproduction, beam-gas 
collisions and cosmic rays is negligible. 

The final data samples consist of (i) 739 events for an integrated luminosity of £ = 
0.418 (±4.5%) pb-i of the e- p 1993 data, (ii) 810 events for C = 0.491 (±2.4%) pb"^ of 
the e-p 1994 data, and (iii) 5201 events for C = 2.947 (±1.8%) pb~^ of the e+ p 1994 data. 

The measured cross sections da/dQ"^ are corrected for detector effects and QED radiation 
and extrapolated to the full kinematic phase space. The correction factors for each bin are 
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obtained from Monte Carlo studies. Deep inelastic scattering events are generated according 
to the Standard Model cross section 
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using the parton densities from the MRS H parametrization ||T0| of the structure functions 
Fi{x,Q'^). The longitudinal structure function Fl has not yet been measured at HERA, but 
is expected to give very small contributions only at low and low x values. For the present 
analysis Fl has been neglected. Radiative effects are taken into account by the event gene- 
rator DJANGO 6 which includes the 0{a) electroweak corrections and the QCD matrix 
elements to first order in a^, supplemented by leading-logarithmic parton showers. The lepton 
detection and measurement is simulated by smearing the event vertex and the generated four- 
vector according to measured resolutions, acceptances and efficiencies, as determined from 
data. The hadron final state is simulated by smearing the four-vector of the 'struck quark' 
as calculated from the electron kinematics (no hadronization) . This very fast and efficient 
acceptance simulation is completely adequate to describe all properties of the deep inelastic 
scattering events used in the present inclusive analysis. 



5 Results 

Cross Sections 

The corrected differential cross sections da/dQ'^ for the e~ p and e'^ p data are shown in Fig. |l] 
and hsted in Table |l[ Statistical and systematic errors are added in quadrature, except for an 
overall normalization uncertainty. The systematic errors include uncertainties of the trigger and 
detection efficiencies, the vertex reconstruction and the lepton energy calibration of the LAr 
calorimeter (all determined from the data), as well as uncertainties due to radiative corrections 
and the choice of different parton distributions (MRS H, MRS D*^' and MRS D~'). The syste- 
matic uncertainties are typically ~ 2%, except for the energy calibration, which contributes 
~ 6% in the low Q"^ region. The overall normalization errors are 3.5% for the e~ p data and 1.8% 
for the p data and account for uncertainties of the luminosity measurement. The acceptance 
is a smooth function of the squared momentum transfer. It rises from ~ 65% at the lowest 
to ~ 80% at ^ 1000 GeV^ and then slowly decreases to ~ 60% at ~ 10, 000 GeV^. The 
extrapolation to the full phase space is not critical for the low and medium region, where 
the validity of the Standard Model is well established. At high Q"^ and large y the Standard 
Model is assumed. The effects of the extrapolation on the contact interaction analysis are, 
however, small. 

The measured cross sections are well described by the Standard Model expectations over 
five orders of magnitude in the range between 160 GeV^ and 20,000 GeV^. No significant 
deviation is observed for either lepton charge, see Fig. 

It is interesting to note that the e~ p cross section tends to be slightly higher than the e"*" p 
cross section at high values, as expected from the different couplings of the leptons to the Z 
boson. The charge asymmetry is, however, not yet significant due to the limited electron data 
statistics. 

The contact interaction analysis investigates the differential cross sections da/dQ"^ and in- 
terpretes any deviation from the Standard Model as lower limits on the ratio mass over coupling 
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Q"^ (GeV^^) (GeV^) 



Figure 1: Differential cross sections dajdQ^ versus (•) in comparison with the Standard 
Model expectations ( — ). a) and b) e^p data, c) and d) e+p data. The error bars repre- 
sent statistical and systematic errors added in quadrature, except for an overall normalization 
uncertainty of 3.5% (e~p data) and 1.8% (e^p data), respectively. 
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Table 1: Differential cross sections da/dQ"^ for the e~ p and p data with statistical and 
systematic errors. There is an additional overall normalization uncertainty of 3.5% {e~ p data) 
and 1.8% (e+p data), respectively, due to systematic errors of the luminosity measurement. 

strength of new heavy leptoquarks or composite e q structures. It is assumed that possible de- 
viations are caused by only one new boson exchange at the time. A combined analysis of the 
e~ p and e"*" p data is performed including for each data set an individual overall normalization 
constant fnorm with its corresponding error Anorm (see discussion above) 

^'-'^k iQ"^) fnorm J \ ^norm J J 

The inner sum extends over the bins of one data set, while the outer sum is taken over both the 
e~ p and e+p data. Thus, fitted parameters are the rjij coefficients of the contact interaction 
model and the two normalization constants. Limits at 95% confidence level are derived from 
the increase of by 3.89 with respect to its minimum value, which in most cases coincides 
with the Standard Model fit. Enlarging the normalization errors arbitrarily by 2% would lower 
the resulting limits by ~ 5%. Choosing different parton distributions, MRS D°' or MRS D~', 
in the cross section calculation changes the limits by ~ 3% in either direction. 

Leptoquarks 

The results of the leptoquark analysis are summarized in Table ^. Only those lower limits on 
Mlq/X are quoted, which exceed the kinematic phase space of HERA for direct production 
assuming a strong coupling of A = 1. 

One notices that vector leptoquarks with a positive coupling to u quarks provide the most 
restrictive bounds approaching O {1 TeV). The sensitivity to scalar leptoquarks is generally 
lower by a factor of ~ 2 and two of them, Sq and S^^2: have only couplings to d quarks. It is not 
obvious which leptoquarks will provide the most stringent limits, when arguing alone on the 
basis of the assigned fermion numbers F and quark densities in the proton. In addition to the s 
channel amplitudes of direct production the contact interaction ansatz implicitely contains the 
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leptoquark coupling to u quark coupling to d quark F Mlq / A 
[GeV-2] [GeV-2] [GeV] 
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cR 

^1/2 

oL 

•^1/2 



Vrr 
VIr 

vIl 



Vrr 



vIr - 



-1 (A/M^q)^ 



2 
2 

2 



2 



350 



360 
340 



^0 








= +1 {X/Mlq)' 
= +1 (X/Mlq)' 







yR 






Vrr 







yR 


Vrr 


= +1 {X/Mlq)' 









760 








viR 


= +1 (X/Mlq)' 


2 


300 




Vrl 


= +1 (X/Mlq)' 


Vrl 


= +1 (A/M^q)2 


2 


710 


^1/2 


VIr 


= +1 {X/M^Qf 
= +2 (X/Mlq)^ 






2 


800 






vIl 


= +1 {X/Mlq)' 





1020 



Table 2: Contact interaction coefficients 77^^, fermion number F and lower limits at 95% con- 
fidence level on Mlq/X for scalar (S) and vector (V) leptoquarks. The leptoquark notation 
indicates the lepton chirality L, R and the weak isospin / = 0, 1/2, 1. The leptoquarks 5" and 
V differ by two units of hypercharge from S and V, respectively. 



crossed t and u channel diagrams. Moreover, not all coupling coefficients rjif contribute with 
the same weight, thus spoiling the naive expectation^. 

As an illustration of the sensitivity of the data to virtual leptoquark exchange Fig. ^ a shows 
the allowed contribution of a vector leptoquark V-^. The contact interaction contributions to 
the Standard Model rise with as expected and amount to ~ 25% at the highest values. 

The indirect limits M^q / X derived from virtual leptoquark exchange nicely complement the 
direct searches [jT], 0, because they exceed the kinematic range of HERA at large couplings. 
The couplings A can be safely extracted down to the nominal centre of mass energy, since the 
present integrated luminosities restrict the accessible masses to M^q = y/xs < 250 GeV. For 
example, a vector leptoquark with a mass of 300 GeV can be excluded for couplings larger 
than the electromagnetic strength A > \^Att a = 0.3. 

The limits given in Table |^ can be compared to those derived from other, primarily low 
energy experiments. The strongest bounds for leptoquarks coupling to the ffist lepton and quark 
generations arise from atomic parity violation experiments and from universality in leptonic tt 
decays. Davidson et al. [|12| give Mlq/X limits in the range of 1.6 to 2.2 TeV for scalar 
leptoquarks and 2.2 to 3.1 TeV for vector leptoquarks, while Leurer |1^ arrives at values up 
to a factor of ~ 1.5 either lower or higher, and generally gets larger values for scalar than 
for vector leptoquarks. Despite of these uncertainties, our most stringent limits for vector 

comprehensive general discussion on the sensitivity of the couphng coefficients in the eight-dimensional 
rji f space can be found in ref . Q . 
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(GeV^^) (GeV^) 

Figure 2: Differential cross sections da/dQ"^ normalized to the Standard Model expectation 
versus for the combined e~ p and p data (•). The error bars represent statistical and 
systematic errors added in quadrature, but do not contain an overall luminosity uncertainty 
(see text), a) The curve shows the allowed contribution (95% confidence level) of a vector 
leptoquark V-t with M^q/A = 1020 GeV. b) The curve shows the effect of a form factor with 
a quark radius Rq = 2.6 • 10^^^ cm (95% confidence level). 



leptoquarks are within a factor of two close to the ones extracted from the very low energy 
experiments. They provide, however, very useful complementary information at much higher 
momentum transfers from a model independent analysis involving no theoretical assumptions 
on higher order corrections. 

Compositeness Scales 

If quarks and leptons have a substructure and have common constituents they may form com- 
posite objects. Such virtual states are characterized by a compositeness scale parameter A and 
a coupling strength g, which is set to (7^/4 tt = 1 in the present analysis. 

The results for lower limits on the eg compositeness scale parameters are summarized 
in Table 0. They vary between 1.0 TeV and 2.5 TeV, depending on the chiral structure and the 
sign of interference with the Standard Model currents. The bounds with positive interference 
are more stringent than those with negative interference. There is almost no difference between 
various lepton and quark chiralities. 



A 



LL 



A 



LR 



A 



LR 



A 



RL 



A 



RL 



A 



RR 



A 



RR 



2.3 1.0 2.5 1.2 2.5 1.2 2.3 1.0 



Table 3: Lower limits at 95% confidence level on compositeness scale parameters A^ [TeV] 
for chiralities LL, LR, RL and RR with positive and negative interference with the Standard 
Model currents. The coupling constants are igf^/A^^ with the convention g^/4:TT = 1. 
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Although it is common practice to give bounds on A^, it is more appropriate to fit directly 
Tjif = g"^ / as used in the contact interaction Lagrangian. This has the advantage that the 
errors behave rather gaussian. Limits on A^*" for a certain interference are then derived from 
the corresponding upper and lower values of 1/A^ at a given confidence level. A fit of 1/A^ 
yields Kll = A^l = (-0.41 ± 0.33) TeV'^ and A^^ = A^l = (-0.30 ± 0.23) TeV'^ assuming 
a coupling strength of (7^/4 7? = 1. 

The compositeness scale parameters A can be interpreted in terms of a radius of the eg 



system via Reg = y4:7i/g''^ A ^. Depending on the chiral structure and the sign of interference 
the size of a composite eg state is constrained to i?eg ^ (0.8 -=-2) ■ 10~^^ cm, if a strong coupling 
strength g' is assumed. 

Our limits on A^ are comparable to those derived in similar studies a.t pp and e"*" e~ col- 



liders IjTj]. For example CDF quotes A'l^ > 1.7 TeV and A^^ > 2.2 TeV for the light u and d 
quarks and VENUS quotes A^j^ > 1.2 TeV and A^^ > 1.6 TeV assuming flavour universality 
for five quarks. 

Form Factors 

An alternative method to study possible fermion substructures is to assign a finite size of 
radius R to the leptons and/or quarks [0. A convenient parametrization is to introduce in a 
'classical' way form factors /(Q^) at the gauge boson-fermion vertices, which depend on the 
squared momentum transfer 



For simplicity, the radius R is assumed to be universal for the electromagnetic and the weak 
vector and axial-vector fermion couplings. The form factors reduce to unity and the couplings 
to the familiar Standard Model values in the pointlike limit. A finite extension of a lepton or 
quark is expected to diminish the Standard Model cross section at high according to 

da da^'^' 



dQ^ dQ^ 

The data are analyzed in terms of a single form factor, yielding as an upper limit at 95% 
confidence level a radius of 

R < 2.6 ■ 10"^^ cm . 

This result may be interpreted as a limit on the light quark sizes, since the pointlike nature 
of the electron is already established down to much lower distances in e"*" e~ and {g — 2)e 
experiments [jT6[. Fig ^b shows the effect of a form factor with a quark radius Rq given by the 



experimental limit on the differential cross sections; again the sensitivity rises with Q^. 

The limit on Rg is within a factor of two comparable to the bounds derived from a global 
analysis of Z decays at LEP [|1^ 



Similar upper limits on a quark radius (up to a factor of v^) are obtained from the above 
contact term analysis, if the compositeness scale parameters A" are evaluated at the electro- 
magnetic scale g'^/'iir = a. Note, however, that in this interpretation the size is inferred from 
the interference of the Standard Model currents with a new virtual current. 



13 



6 Conclusions 



An analysis of searches for new phenomena beyond the Standard Model mediated through con- 
tact interactions in deep inelastic scattering at HERA has been presented. The data correspond 
to a 9 fold increase in integrated luminosity compared to a previous publication [|l| . 

The differential cross sections da jdQ^ have been measured for deep inelastic neutral current 
e^p scattering in the range between 160 GeV^ and 20,000 GeV^. No significant deviations 
from the Standard Model have been observed for either lepton charge. 

Substantially improved limits at 95% confidence level on masses and couplings of new heavy 
leptoquarks and on fermion compositeness scales have been obtained, using the combined e~ p 
and p cross section data. 

Eight out of fourteen possible leptoquark couplings yield lower limits on Mlq/X which 
exceed the centre of mass energy of HERA, assuming a strong coupling A = 1, and therefore 
nicely complement the searches for direct production. Vector leptoquarks yield stronger limits 
than scalar leptoquarks and approach bounds on M^q/X of 1 TeV. 

A conceivable eg compositeness or fermion substructure can be ruled out for scale para- 
meters A^ smaller than 1.0 TeV to 2.5 TeV, depending on the assumed fermion chiralities and 
the sign of interference with the Standard Model currents. 

Finally, a form factor analysis constrains the size of the light u and d quarks to radii smaller 
than Rg < 2.6 ■ lO-^^ cm. 

The contact interaction concept has been shown to be a very powerful tool, becoming even 
more important in future high statistics data analyses. The sensitivity to new virtual boson 
exchanges roughly scales as {C s)^, with integrated luminosity C and centre of mass energy 
squared s. 
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